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Abstract - The Letter addressees vortex dynamics for the lattice models of superfluids in the 
continuous approximation. The continuous approximation restores translational invariance absent 
in the original lattice model, but the theory is not Galilean invariant. As a result, calculation of 
the two transverse forces on the vortex, Magnus force and Lorentz force, requires the analysis of 
two balances, for the genuine momentum of particles in the lattice (Magnus force) and for the 
quasimomentum (Lorentz force) similar to that in the Bloch theory of particles in the periodic 
potential. While the developed theory yields the same Lorentz force, which was well known before, 
a new general expression for the Magnus force was obtained. The theory demonstrates how a small 
Magnus force emerges in the Josephson-junction array if the particle-hole symmetry is broken. The 
continuous approximation for the Bose-Hubbard model close to the superfluid-insulator transition 
was developed, which was used for calculation of the Magnus force. The theory shows that there 
is an area in the phase diagram for the Bose-Hubbard model, where the Magnus force has an 
inverse sign with respect to that which is expected from the sign of velocity circulation. 



Introduction. — The transverse force on a vortex in 
superfluids (neutral and charged) is debated during many 
decades and has been a topic of reviews and books [lr|4|. 
In a continuous supcrfluid at T = 0, which is identical to 
a perfect fluid in classical hydrodynamics, the balance of 
forces on a vortex is 



F L = mn \{v L - v s ) x k] = F ext , 



(1) 



where the Magnus force Fm is proportional to the vor- 
tex velocity Vl and the Lorentz force Fl is proportional 
to the superfiuid velocity v s = ^V<^ determined by the 
phase tp of the order parameter wave function, and the 
external force F ext combines all other forces on the vor- 
tex, e.g., pinning and friction forces. Here n is the density 
of particles with mass m, and k is the vector parallel to 
the vortex axis with its modulus equal to the circulation 
quantum k = h/m. The united transverse force -Fm + Fl 
depends only on the relative velocity vl — v s as required 
by the Galilean invariance. 

In lattice models of superfluids the Galilean invariance 
is absent, and the value of the Magnus force was under 
scrutiny. The most known lattice model of the superfiuid 
is the Josephson junction array. Usually they studied the 
theory of vortex motion in the array using the continu- 



ous approximation. These studies have not revealed any 
Magnus force normal to the vortex velocity j5]. Moreover, 
there have been experimental evidences of the ballistic vor- 
tex motion in the Josephson junction array [5], which is 
possible only in the absence of the Magnus force. There 
was a lot of theoretical works trying to explain this impor- 
tant feature of the Josephson junction array, mostly sug- 
gesting that quantum effects can provide a finite Magnus 
force. Meanwhile, in the classical theory of the Josephson 
junction array they usually assumed the particle-hole sym- 
metry, which forbids the Hall effect and the Magnus force 
in the model [7]. However, this symmetry is not exact, and 
its violation allows the existence of a finite Magnus force. 
In superconductors the Magnus force determines the Hall 
effect, and the presence or the absence of this force means 
the presence or the absence of the Hall effect. 

Intensive experimental and theoretical investigations of 
Bosc-condcnsed cold atoms attracted an interest to an- 
other lattice model of a superfiuid: the Bose-Hubbard 
model [8]. The periodic structure of potential wells for 
bosons, which leads to the Bose-Hubbard model in the 
tight-binding limit, is realised for cold-atom BEC in the 
experiments with optical lattices [9]. Recently Lindner et 
al. [TU] and Huber and Lindner [TT] addressed the prob- 
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lem of the Magnus force in the Bose-Hubbard model and 
revealed that close to the superfluid-insulator transition 
the force changes its sign as happens in Fermi supcrfluids 
when the sign of the carrier charge changes its sign. 

The Letter presents the analysis of the transverse force 
on the vortex in a lattice, which is approximated by a 
continuous model with effective parameters determined by 
properties of the lattice. The forces are determined from 
the momentum balance. The absence of the Galilean in- 
variance makes necessary to analyse two momentum bal- 
ances: for genuine momentum and for quasimomentum 
similar to that in the Bloch band theory for particles in 
a periodic potential. This yields the general expression 
for the Magnus force in the absence of the Galilean in- 
variance. The expression allowed to calculate the Magnus 
force in the Bose-Hubbard model close to the superfluid- 
insulator transition after deriving the continuous theory 
for the supcrfluid order parameter near this transition. 

Vortex dynamics in the continuous approxima- 
tion for the lattice superfluid. — The continuous ap- 
proximation for lattice superfluids gives the phenomeno- 
logical theory reminding the usual hydrodynamics. Gen- 
erally enough this approximation corresponds to the fol- 
lowing Lagrangian: 



L 



-hntp 



h 2 n 
2m 



(Vcp) 2 - E c {n). 



(2) 



where E c (n) is the energy of a resting liquid which de- 
pends only on n. For simplicity we consider the 2D prob- 
lem, where n is the particle number per unit area. The 
Hamiltonian (energy) for this Lagrangian is 



H 



dL 



dip 



L 



Ti 2 n 
2m 



(V<p) 2 + E c (n). 



(3) 



Despite similarity to hydrodynamics of the perfect fluid, 
there is an essential difference. The continuous approxi- 
mation for the lattice model restores translational invari- 
ance but not Galilean invariancc. The latter is absent 
since the effective density h, which characterises stiffness 
of the phase field, is different from the genuine particle 
density n, and is much less than n if the lattice nodes are 
weakly connected. 

Let us discuss the conservation laws, which follow from 
Noether's theorem. The gauge invariance provides the 
conservation law for charge (particle number): 



d dL „ ( 8L 
h V fe 

dt dip \9Vfe<y9 



= 0, 



(4) 



which in fact is the continuity equation (the first Hamilton 
equation) for the density n: 

dn h _ .__ , ._. 

-77 = V fe (nV^). (5) 

dt m 

The second Hamilton equation is 

dip dH H 2 dn 2 

n m = -dn- = -»-2^^ j(p) ' (6) 



where u = dE c (n) / dn is the chemical potential of the 
liquid at rest. 

The translational invariance provides the conservation 
law 



d fdL^ \ _ ( dL 

dl\M v ^) +Vl [dV^ 

This is the conservation law 
dgk 



\7ip>- L5 t 



0. 



dt 



+ v ; n fe/ = 0, 



for the current 



g = -—Vip = hnvtp. 
Here the momentum-flux tensor is 



m 



(7) 



(8) 



(9) 



P+ ^[tn n -^^^ 



'ki- 



rn 



and the pressure is determined from the T = thermody- 
namic Gibbs-Duhem relation dP = ndn. 

The current g is not a genuine mass current. The gen- 
uine mass current is that appears in the continuity equa- 
tion ([5]) and is the density of the genuine momentum: 



3 = 



dL 



h dVp> 



= HnVip. 



(11) 



So Noether's theorem does not provides the conservation 
law for the momentum in a non-Galilean liquid. Neverthe- 
less, the momentum conservation law for the Lagrangian 
@ takes place though it is derived not from Noether's 
theorem but directly from explicit equations of motion: 



djk 
dt 



0, 



where the momentum-flux tensor is 
h 2 dh 

n fc i = ——nVk^iip + PSki, 
m an 



(12) 



(13) 



and the partial pressure P is determined from the relation 
dP = hdu. These equations reduce to those in hydro- 
dynamics of a Galilean invariant liquid when h = n. In 
contrast to the conservation law eq. © following from 
Noether's theorem, the conservation law (TTZj) for the gen- 
uine mass current is approximate and neglects higher than 
second order terms in phase gradients. 

For better understanding of the physical difference be- 
tween the two currents g and j, let us consider particles 
in a periodic potential, which can be also described by 
our continuous model. The wave functions of particles are 
described by Bloch functions: 



ip(r) = u n (r,k)e 



ih-r 



(14) 
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where u n {r, k) is the periodic function and the wave vector 
k determines the quasimomentum hk. The quasimomen- 
tum differs from the genuine momentum of the quantum 
state equal to the average of the momentum operator: 



Then the momentum-flux tensor (jlOp becomes 



p = —ih 



hk — ih 



or 



(15) 

Calculating the band energy E(k) in the kp approxima- 
tion for small k, i.e., at the band bottom (the energy min- 
imum at k = 0) one obtains that 



P 



where 



dE{k) 



-hk, 



d 2 E(k) , hk 
-k = , 

m 



(16) 



(17) 



9 hdk hdk 2 
is the group velocity and m* is the effective mass. The vec- 
tor p is a genuine momentum indeed since particles with 
bare mass m propagate in space with the group velocity. 
Assuming that particles are bosons, which condense in a 
single Bloch state with density n, and connecting the wave 
vector k = V<p with the phase tp, the density of the gen- 
uine momentum, i.e. the genuine mass current is, j = np 
coincide with that given by eq. (|TTj) with h = nm/m*. On 
the other hand, the current g = nhk is the quasimomen- 
tum density. It is well known from the solid state physics 
[12] that an external force on a particle in the energy band 
determines time variation of the quasimomentum but not 
the genuine momentum: 



dk t dv g 
h— = m —r- = f. 

dt dt 



(18) 



In the absence of Umklapp processes the total quasimo- 
mentum is also a conserved quantity, and the conservation 
law for the Bose-condensate in a single Bloch state is given 
by eq. ([8]). The physical meaning of the effective second 
Newton law (fT8| is that only the part dp/dt of the whole 
momentum hdk / dt brought to the system by the external 
force is transferred to to the momentum of particles mov- 
ing in the periodic potential. The rest part d(hk — p)/dt 
is transferred to the system supporting the periodic po- 
tential, i.e., to the lattice. 

For derivation of forces on the vortex one must calcu- 
late the momentum flux through a cylindric surface sur- 
rounding the vortex line [TJE] ■ We argue that the Lorentz 
and the Magnus force must be derived from the balance 
of different momenta: the quasimomentum for the former 
and the genuine momentum for the latter. Deriving the 
Lorentz force one can assume that the vortex is at rest in 
the laboratory system connected with the lattice. Solv- 
ing eq. ([6]) for the time-independent phase tp one obtains 
the quadratic in Vy? correction to the chemical potential 
(Bernoulli's effect): 



h 2 dh ,_ , 5 
2m dn 



(19) 



n» 



m 



[Po - n(V,^) 2 ] 4; 



(20) 



where Po is a constant pressure in the absence of any 
velocity field. The components of the Lorentz force are 
given by the integral over the cylinder around the vor- 
tex: F^i = ifllkidSi. Taking into account that the 
phase gradient Vy> = Vp v + Vtpj consists of the gradient 
Vi^u = [z X r]/2irr 2 induced by the vortex line and the 
gradient 'Vipj = j/hn produced by the transport current, 
the integration yields the Lorentz force 



-[j x re] 



-mn u, x re 



(21) 



This expression for the Lorentz force is well known but 
now we derived it from the quasimomentum balance. This 
is because the force is a transfer of momentum from the 
transport velocity field to the vortex. But any variation 
of the transport velocity must be accompanied by the mo- 
mentum transfer to or from the lattice. So as in the effec- 
tive second Newton law (fT8|) . the force is a time variation 
of quasimomentum but not of genuine momentum. 

The derivation of the Magnus force requires the balance 
of the genuine momentum. It is more convenient to con- 
sider this balance in the coordinate frame moving with 
the vortex where variables are time-independent. In the 
frame moving with the velocity w the energy E' is con- 
nected with the energy E in the laboratory frame by the 
relation E' = E — j ■ w (we ignore quadratic in w terms). 
Then according to the Hamiltonian ([3]) written in the lab- 
oratory frame 



E' 



h n , 
2m 



(Vip') 2 + E c (n), 



(22) 



where the phase gradient in the moving frame is Vy/ = 
Vip—mw/h. This does not differ from the Hamiltonian in 
the laboratory frame, and one can use the momentum-flux 
tensor (|13[) replacing tp by ip'. The pressure variation is 
determined from the Bernoulli law as before. Deriving the 
Magnus force one may ignore the superfluid current in the 
laboratory frame, i.e. Vyj = and replace w by v^. Then 
calculating the Magnus force components F^ji = § ^-udSi 
one obtains the Magnus force 



_ dn _ . 
F M = —nm[v L x re . 
dn 



(23) 



The force appears due to convection of the vortex-related 
momentum hnV<p v into the area of the momentum bal- 
ance by the supercurrent. Since the circular velocity field 
^ Vtp v is fixed an arrival of a particle into the balance area 
may change the vortex-related momentum only via varia- 
tion of n and does not require accompanying momentum 
transfer to the lattice. This is why the Magnus force is 
determined from the balance of the genuine momentum 
and is proportional to dh/dn. 
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In the Josephson junction array the current between two 
nodes of the lattice is determined by the Josephson energy 
Ej cos(<pi — (^2), where tpi and if2 are the phases of the two 
nodes. In the continuous limit this yields n = mEj/K 2 . 
The particle-hole symmetry requires that Ej and n do not 
depend on the average density n, and the Magnus force 
vanishes in agreement with the symmetry of this model. 

Including the Lorentz and the Magnus force into the 
balance of forces on the vortex, like that given by eq. ([T]) 
for the Galilean invariant liquid, one see that the total 
transverse force Fm + Fl is proportional to the velocity 
(dh/dn)vL — v s , but not to the relative velocity vl — v s . 

The role of the first term —hntp in the Lagrangian 
which is called the Wess-Zumino term, was widely dis- 
cussed in the literature [4]. The arguments were about 
whether the total liquid density n must be replaced by 
some other density. It is evident that adding any con- 
stant C to the density n in the Wess-Zumino term docs 
not affect the dynamical equations. However the role of 
the Wess-Zumino term changes after transition from the 
continuous model in terms of fields to the reduced descrip- 
tion in terms of the vortex coordinates r^{xL, Vl)- For the 
Wess-Zumino term this leads to substitution of the phase 
field ip v {r,ri,) for moving vortex into the term and its 
integration over the whole space. Bearing in mind that 
ip v = —(vl ■ V)(/3„, the Wess-Zumino term becomes 



Lwz = -h{n + C)v L ■ [z x r L }. 



(24) 



Varying the total Lagrangian of the vortex with respect 
to r^it), one obtains the equation of vortex motion with 
the effective Magnus force oc (n + C). So the constant 
C does matter for the value of the Magnus force. Some- 
times it is argued that the contribution oc C is of topolog- 
ical origin and can be found from the topological analysis 
[J. On the other hand, we think that there is no general 
principle, which dictates the charge in the Wess-Zumino 
term. Not the undefined Wess-Zumino term determines 
what and whether the transverse force is, but vice versa, 
one must derive the transverse force from dynamical equa- 
tions and only after this one knows what Wess-Zumino 
term in the vortex Lagrangian should be. In particular, in 
the Galilean translation invariant liquid one obtains from 
the momentum-conservation law that the amplitude of the 
Magnus force is proportional to the total density. Then 
only the latter enters the Wess-Zumino term and C = 0. 
On the other hand, if in the Fermi supcrfluid the core 
bound states are intensively scattered by impurities, the 
Kopnin-Kravtsov force [5] fully compensates the Magnus 
force. Then the vortex Lagrangian does not contain the 
Wess-Zumino term at all. 

Vortex dynamics in the Bose Hubbard model. 

— The Hamiltonian of the Bose-Hubbard model [5] for a 
lattice with distance a between nodes is 

h = -jY. s & + j E Kin -i)-mE^- ( 25 ) 



Here /j, is the chemical potential, the operators hi and h\ 
are the operators of annihilation and creation of a bo- 
son at the ith lattice node, and Ni = b\bi is the particle 
number operator at the same node. The first sum is over 
neighbouring lattice nodes i and j. 

In the supcrfluid phase with large numbers of particles 
Ni all operator fields can be replaced by the classical fields 
in the spirit of the Bogolyubov theory: 



'Nte 



(26) 



where <pi is the phase at the ith node. After transition to 
the continuous approach one obtains the Hamiltonian ([3]) , 
whcrc0 



N 



mzoJa' 



n, E c (n) 



Ua 2 



fxn. (27) 



Here zo is the number of nearest neighbours equal to 4 
in the quadratic lattice. Thus the effective density n is 
proportional to the total density n, as in the case of the 
condensate of bosons in a single Bloch states close to the 
bottom of the energy band considered in the previous sec- 
tion. This is not accidental since in the limit of weak in- 
teraction U the hopping term ~ J in the Hamiltonian also 
yields the energy band, which is described by the effective 
mass m* near its bottom, and the ratio h/n is equal to 
the ratio of masses m/m* . According to derivation of the 
previous section, the Magnus force [see eq. (|23p ] is by the 
factor n 2 /n 2 = m 2 J 2 a' l /h i smaller than in the Galilean- 
invariant supernuid. This factor must be small because of 
weak Josephson links between the potential wells. 

It is known that when the energy J of the internodc 
hopping decreases, the phase transition from supcrfluid 
to Mott insulator must occur [8]. In the limit zqJ/U —¥ 
when the hopping term oc J can be ignored the cigenstates 
are given by Fock states |$jv) = \N) with fixed number N 
of particles at any node. At growing J the transition line 
can be found in the mean- field approximation |13) . In this 
approximation the hopping term is taking into account by 
introducing the mean field equal to the average value of 
the annihilation operator its complex conjugate being the 
average value of the creation operator: 



<6. i )=V^ = |^|e^, (b\) = i>t = Me"**" 



(28) 



It is assumed that only the phase but not the modulus 
of the order parameter tp varies from node to node. In 
contrast to eq. ([26)) . in general \tp\ 2 is not equal to N 
and must be determined from the the condition of self- 
consistency (sec below). Introducing the mean field one 
reduces the problem to the single-node problem with the 
Hamiltonian 



U a = -zJ(b\ip + + ^Ni(Ni - 1) - uNi. (29) 



U 



1 We assume that there is the same number of particles at all 
nodes and write Ni without the subscript i. 
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Superfluid phase 



Region of the inverse 
Magnus force 



TV = 



z J/U 



Fig. 1: The phase diagram of the Bose-Hubbard model. The 
Mott insulator phase occupies lobes corresponding to fixed in- 
teger numbers TV of bosons. The shaded beaks of the superfluid 
phase, which penetrate between insulator lobes, are analysed 
in the text. The dash line separates the region with the inverse 
Magnus force from the rest of the superfluid phase. The line is 
schematic since it was not really calculated. The region of the 
inverse Magnus force exists under any lobe but is shown only 
for the beak between the N = 1 and N — 2 lobes. 



Here 



E' 



(30) 



reduces to the number zq of nearest neighbours in the 
uniform state with the constant phase at all nodes. 

The many-node wave function is a product of the single- 
node wave functions. Calculating the energy of the origi- 
nal Hamiltonian (|25[) for this wave function and minimis- 
ing it with respect to ip one obtains the self-consistency 
equation, which determines ip. Following this approach 
|13j one obtains the phase diagram shown in Fig. [TJ The 
Mott-insulator phases with fixed numbers N of particles 
per node occupy the interiors of lobes at small zqJ/U. 

We are interested to know what vortex dynamics is ex- 
pected close to the phase transition at minimal values of 
J, i.e. at beaks of the superfluid phase between lobes, 
which are shaded in Fig. [TJ Here the mean-field approxi- 
mation is simplified by the fact that only two states with 
N and N + 1 particles interplay at the superfluid phase in 
the beak corresponding to the given N. This is because at 
\i = NU these two states have the same energy, whereas 
all other states are separated by the gap on the order of 
the high energy U. So for a beak between the lobes cor- 
responding to Mott insulators with the number of bosons 
per node N and N + 1 we look a solution in the form of a 
superposition of two Fock states: 



\^ N ) = fN\N) + f N+1 \N + l). 



(31) 



This wave function is an eigenfunction of the single-node 
Hamiltonian (|29p if 



TN 



lf-+Z 2 J 2 {N +1 



2 _ H_ 
2 



In+i = 



U£ + z 2 P(N + 1)H 



2i.£ 
2 



The energy of this cigenstate is 



(32) 



e N = -ft' [ N + - | 



z 2 .P{N + 1)|V| 2 , (33) 



where fi' = fi — Un, and one should choose the lower sign 
- for the ground state. The average number of particles in 
the node is a function of //: 



IN) =N + - + — /jt 

2 2^/i' 2 +4z 2 ,P(n + l)\iP\ 2 



(34) 



The self-consistency equation follows cither from the min- 
imisation of the total energy with respect to ip or from the 
condition that ip is the average value of the operator a: 



= (b) = 



zJ(N+l) 



2 \/4 +z2j2 ( iV+1 )IVf 



-.ip. 



(35) 



A non-trivial (i.e., non-zero) solution of this equation is 



7V+1 



fi 



Az 2 J 2 (N + 1)' 



(36) 



The eigenvalue e^v determines the average value of the 
Hamiltonian (|25l) . which is the Gibbs thermodynamic po- 
tential of the grand canonical ensemble per one node: 



Gn = zJ\ip\ 2 + e N . 



(37) 



Since the Hamilton equations are formulated for the 
canonically conjugated variables "phase-particle number" 
it is useful to go from the grand canonical ensemble with 
the Gibbs potential being a function of \i to the canoni- 
cal ensemble where the energy density is a function of the 
particle number density n. Then the energy per node is 

UN 2 

E N = G N + fiN = — — + UNN e 



+zJ (\ip\ 2 - 2VNTT^-N 2 \iJj 



where 



N e = (N)-N 
The energy has a minimum at 

h/f = (iv + i) / 



n: 



(38) 



(39) 



(40) 
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As is expected in the case of the second-order phase tran- 
sition, ip vanishes at the phase transition lines, where 
A e = ±i and the number of particles reaches A at the 
lower border and N + 1 at the upper one. In contrast 
to the Landau-Lifshitz theory of the second order transi- 
tions, there is no analytic expansion in ip near the critical 
temperature because of the term linear in \ip\. 

In order to derive parameters of the continuous model 
discussed in the previous section, let us consider the effect 
of slow phase variation from node to node. Assuming that 
ipi = \ip\e lk ' Vi where r, is the position vector of the ith 
node, one obtains for the square lattice with the number 
zq — 4 of nearest neighbours: 

z = 2cos(k x a) + 2cos(fc y a) « 4 — k 2 a 2 . (41) 

Bearing in mind the correspondence of the wave vector k 
to the gradient operator V in the configurational space one 
obtains in the continuum limit for small k the Hamiltonian 
Q with 

n=^J(N + l) 

E c (n)=4J(N+l)(±-n 2 e aA, (42) 

where the effective density n e = N e /a 2 differs from the 
total density n = (A + |) /a 2 + n e by a constant and 
constant terms in the energy were ignored. This allows to 
find the density dependent factor in the expression (1231) 
for the Magnus force: 

d £n = jV(A + l) 2 n e Q - n 2 e a^ . (43) 

A remarkable feature of the Magnus force in the beaks of 
the supcrfluid phase is that its sign is not a sign of the 
total charge proportional to the total density of the liquid 
but is a sign of the effective density n e . So in the upper 
halves of the beaks shown in Fig. [T]the Magnus force has 
an inverse sign with respect to the total charge. If the 
liquid were electrically charged this would mean an inverse 
sign of the Hall conductivity. The regions of the inverse 
Magnus force neighbour any insulator lobe from below, 
where n e is positive. Since at upper borders of the lobes 
n e is negative, the line n e = must end somewhere at the 
border of the lobe. In Fig. Q] it is shown by a dashed line 
for the beak between the N = 1 and N = 2 lobes. The 
connection of the Hall conductivity with the sign of n e 
was pointed out by Lindner et al. [10] for the lowest beak 
A = 0, i.e., for the superposition of the states A = and 
A = 1 at the node, which characterises the model of hard- 
core bosons (the limit U — > oo). However they obtained 
the Magnus force essentially different from eq. l|2"3"]) and 
having a jump at the line where the force changes sign. 

Conclusion. — We derived the transverse (Magnus 
and Lorentz) forces on the vortex from the conservation 



laws in the continuous approximation for lattice models 
of superfluids. It is important that the two forces are ob- 
tained from two different conservation laws, one for the 
genuine momentum of particles in the lattice (Magnus 
force), another for the quasimomentum (Lorentz force) 
similar to that in the Bloch band theory The calculated 
Magnus force vanishes for the Josephson junction array as 
required by the particle-hole symmetry usually assumed in 
the theory of the array. The analysis of the Bose-Hubbard 
model in the continuous approximation yielded the ex- 
pression for the Magnus force, which had an inverse sign 
at some area close to the phase superfluid-insulator tran- 
sition. Sign inversion was revealed earlier for hard-core 
bosons by a different method [TOlfTT] but from the expres- 
sion for the Magnus force different from ours. 

Though our analysis addressed ideal strictly periodical 
lattices, the results of the analysis are expected to be valid 
also for superfluids in random potentials (e.g., superfluids 
in porous media or on disordered substates), where broken 
translational and Galilean invariance must also weaken the 
Magnus force. 

* * * 
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